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Employing the self-consistent Green's function approach, we studied the temperature dependence 
of the spin-wave stiffness in diluted magnetic semiconductors. Note that the Green's function 
approach includes the spatial and temperature fluctuations simultaneously which was not possible 
within conventional Weiss mean-field theory. It is rather interesting that we found the stiffness 
becomes dramatically softened as the critical temperature is approached, which seems to explain 
the mysterious sharp drop of magnetization curves in samples within diffusive regime. 
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Spintronics0 brings out an industrial renaissance in 
the past decades for its powerful usage of the extra spin 
degrees of freedom in many electronic applications. For 
instances, giant magneto-resistance (GMR) had been ap- 
plied to read-write head of computer hard disk and tun- 
neling magneto-resistance (TMR) to the newest non- 
volatile memory for MRAM. The huge success kicked 
off the intense investigations on possible realizations of 
similar devices in semiconducting materials, which can 
be directly integrated with the existing industrial tech- 
niques. The idea has enjoyed its primary success in the 
so-called diluted magnetic semiconductors (DMS), where 
the magnetic ions are doped into the host semiconductors 
and lead to a ferromagnetic phase. 

The key issue at current stage is how to raise the crit- 
ical temperature so that the ferromagnetic order is ro- 
bust even at room temperature. Taking the well-known 
material (Ga,Mn)As as exa mpled H H, it was demon- 
strated that the critical temperature can be raised sig- 
nificantly by thermal annealing. However, the highest 
critical temperature at the time of writing is around 160 
Kp|, which is still far from the goal for room-temperature 
DMS. Based on Zener model, Dietl et al. proposed to 
look for room-temperature DMS in wide bandgap semi- 
conductors and oxides such as GaN, ZnO, Ti02 |E @ • 
While the critical temperatures in these materials are 
typically higher compared with (Ga,Mn)As, clustering 
seems to be a serious problem which prevents their po- 
tential usage in realistic devices. 

While the issue of raising up the critical temperature 
seems to lie in experimentalists' hands, we believe a bet- 
ter understanding of the ferromagnetic phase would also 
help. After intense theoretical investigations, the ori- 
gin of ferromagnetism in DMS is believed to be carrier 
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mediated^. Typically, the doped transition metal ions 
provide the localized impurity spins with small direct 
exchange (often antiferromagnetic) among themselves. 
Furthermore, the experimental results show that only 
10^30% of the doped magnetic ions contribute itinerant 
carriers [lfjj into the host semiconducting bands. Through 
double-exchange mechanism, the itinerant carriers medi- 
ated the indirect exchange interactions between the local- 
ized impurity spins and the ferromagnetic order sets in 
when the system is cooled below the critical temperature. 

While the origin of the ferromagnetic phase is more or 
less clear, some of its physical properties remain puzzling. 
To achieve a full understanding, it is crucial to include 
the exchange coupling between itinerant and localized 
spin densities, thermal fluctuations, the random locations 
of the doped ions, the realistic band structure and the 
repulsive interaction between itinerant carriers. Since it 
is almost impossible to incorporate all effects in single 
formalism, one needs to glue up piecewise information 
from different approaches. 

In this paper we employ the self-consistent Green's 
function methodpr|. which includes the spatial and ther- 
mal fluctuations simultaneously, to study the temper- 
ature dependence of spin- wave stiffness in DMS. Since 
our goal is to demonstrate the peculiar temperature de- 
pendence, we simplify the realistic band structure by 
the single-band approximation, at price of sacrificing the 
quantitative description for realistic materials. However, 
within single-band approximation, we were able to estab- 
lish the close tie between the softening of stiffness and 
the sudden drop in magnetization curves^- Moreover, 
our self-consistent Green's function approach also shows 
the appearance of concave magnetization curve |12j in the 
regime where itinerant carriers arc dilute. This implies 
a smooth crossover from the diffusive regime to the lo- 
calized regime with strong disorder. Since previous stud- 
ies including the full six bands only renormalize various 
physical parameters and twist the phase boundaries [l^. 
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we expect our results to be qualitatively robust. 

The organization of the paper is the following: in the 
Section II of this paper, we derive the formalism for the 
self-consistent Green's function approach. In Section III, 
we show our numerical results and discuss connections to 
other approaches in the literature. 

The sd model is proper to describe the DMS systems 
including a strong exchange interaction between local 
spins, which come from electrons in d orbits of transition 
atom, and itinerant spins around whole system, which 
come from impurity doping donation. The Hamiltonian 
of the sd model can be expressed by 



H = H + J / d 3 rS(r) ■ a(r), 



(1) 



the first term Hq is the kinetic energy of the itinerant 
carriers and the second term is the exchange interaction 
between itinerant carrier spins and the localized spin mo- 
ments, where the spin density of the localized moments 
is S(r) — ^2 1 ^ 3 ( r — Rj)Si an d the itinerant spin density 
is er(r) = tjj'(r)(T/2)i/j(r). In momentum representation 
the Hamiltonian is expressed as 



H — Hk + Hi 



Ck,a 



J 



~E 5 ^ + - cJ E 5 ^-*> 



(2) 



where the constant c is the density of the magnetic ions 
in DMS and J is the magnetic coupling integral in the 
unit of eVnm 3 . Since the RKKY interaction dominates 
the magnetism in DMS, to investigate the magnon dis- 
persion is necessary for studying in microscopic. In order 
to obtain the magnon dispersion of local spin in DMS we 
have to calculate the retarded local spin Green's function 
defined by 



Gi,M = ((St(t);Sr(0)y)) 

= -ie(t)([S+(t),Sr(0)}), 



(3) 



where 9(t) is a step function of time t, (• • ■) representing 
the expectation value and [• • •] is the commutator. It is 
more convenient to derive this Green's function in mo- 
mentum space. Through the Fourier transformation to 
obtain the spin Green's function in momentum represen- 
tation, we get 



G(q,t) = ((S+(t);S- (0,0)}). 



(4) 



Employing the equation of motion to equation Q 
makes the spin Green's function calculation reduced to 

ijG{q,t) =<p+({[S+(t),H];S- (0,0))), (5) 

where function ip comes from the derivative of time 
for step function. Following commutation rules in 



RPA (random phase approximation) are we used in this 
paper, there are 



= 2{S z )5 qt k, 



(6) 



(7) 



'p-k M 



and 



(9) 



where (S z ) is the expectation value of magnetization for 
local spin. Employing these commutation rules to equa- 
tion (JSJ results to the equation 

uG(q,w) = <p-cJ(S'){{a+;S-(p,0))) 

+cJ( ( t 2 )«S+;S-(0,0))), (10) 

where (er z ) is the expectation value of the magnetiza- 
tion for itinerant carriers. We found the derivation pro- 
cess for local spin Green's function resulting to another 
new Green's function ({k,uj) = S~(0, 0))) in the 

meanwhile. It is clearly the new resulted Green's func- 
tion reveals the physics that the exchange interaction be- 
tween local spins needs the itinerant carriers' mediation. 
Therefore we need to calculate the new Green's function 
by the same way, 



d 

~dt' 



^<«;S-(0,0)» = ({[a+, H}; S~ (0,0))). 



The commutator in equation l|ll|) results to 

WqiHk] = E( £ P ~ £ P+<l) C l+q,\ C Pd 



(11) 



(12) 



k k 

(13) 

where Sk is the electrons kinetic energy with momentum 
k. Apply the Fourier transformation and RPA for carriers 
density to the new Green's function £ resulting to 



wC(g,w) = E( £ p _£ p+9)(( c p+9,t c p^ ;5 (°'°))) 



C.J 



xG(q,L0)+cJ(S z )C(q,L0). 



(14) 
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Abstracting momentum p from ((q, u>) in equation 1|14|) 
results to a relation 



H+^dlS- (0,0))) 



xG(q, lo). 



cJ(S z ) 
(15) 



Combining equations lfTU)l and (fH|l results to a closed 
form of Green's function equation G(q,u), 



-P+9 



cJ{S* 



xG(q,uj) = tp, 



(16) 



where the (cv a c v ^) — f Pt<7 — {(3e p . a + 1) _1 is the carrier 
density for spin a and /? = 1/KbT. The poles in eq.JTBJl 
represent magnon excitations. In the diluted limitation, 
the minority is the itinerant carriers which could be seen 
as a free carrier gas with effective mass m* — O.5m e 0, 
where m e is the free electron mass, and in the ferromag- 
netic state the majority of local magnetic ions produce an 
effective magnetic field causing both spins to split with 
Zeeman energy cJ{S z ). The kinetic energies of different 
spins are e^^ — h 2 k 2 /2m* =F crl/2cJ(S z ) respectively. 
Therefore the magnon excitation energy for each momen- 
tum q is 



km*fkf 



ln( 



cJ(S z ) 





2-2 



2m 



fc2 

n q 

h 2 kq 



- cJ(S z ) 
km*f kl 



h 2 q 2 
2m* 



h 2 kq 
m* 



)dk + 



*2 

n q 



ln(- 



cj(s z ) + ^4 + r ^' 

\ ' 2m* m* 



(17) 



Finally, we could utilize the Callen's arbitrary spin for- 
mula, 



(S z ) = 



[S - $(S)] [1 + $(g)] 2S+1 + [5 + 1 + $(5)] : 

[1 + $(5)] 2S + 1 - [$(^)] 2S +! 
MS)] 2S+1 



(18) 



to obtain the magnetization values, where $(5) = 
jj £^ g (exp(/3cJ g ) — 1) _1 is the magnon number. 

The anomalous temperature dependance of magnetiza- 
tion of sharp drop in the vicinity of Tc represents in our 
theoretical result as shown in the insert of Fig. 1. is con- 
sistent with experimental result implying an important 
interaction existing in the magnetization collapse region. 
According to the RKKY mechanism the magnetism is es- 
tablished by itinerant carriers mediation. The sharp drop 
of magnetization reveals a possibility of carrier-magnon 
decoupling in the vicinity of Tc- This decoupling effect 
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FIG. 1: We take the exchange coupling and the effective mass 
are fixed at typical values J = 0.15 eV nm 3 , m* = 0.5m e and 
the ratio of itinerant and localized spin densities fixed at c*/c 
=0.1 to calculation resulting to Tc = 45.K". In the insert of 
the figure shows a sharp drop of magnetization in the vicinity 
of Tc- In the sharp drop region the magnon shows a softening 
effect. 



reduces the effective magnetic interaction and the mag- 
netization disappears in decoupling completely. 

The dispersion of spin wave derived from the conven- 
tional spin wave theories|l4j is temperature independent, 
which is an intrinsic characteristic for many kinds of mag- 
nets. The Fig. 1 exhibits the dispersion of magnon 
from our theoretical calculation showing a temperature 
independent dispersion at far from Tqs, which reveals 
a result hat the temperature independence of normal 
spin wave existing at robust magnetism region, mean- 
while it reveals an obvious magnon softening effect in the 
vicinity of Tc leading to the magnetization falling down 
sharply. Interestingly, this softening effect starts from 
small magnon momentum qs then extending to whole 
dispersion region eventually. From the conventional spin 
wave theory as the q -C, the magnon dispersion relation 
has u>(q) = Dq 2 , where the stiffness constant D oc J and 
J is the magnetic coupling integral between two sepa- 
rated spins. From the linear response theory we have 
derived 15] before, the coupling J is oc J 2 , where the J 
is the coupling between itinerant spins and local spins. 
Therefore this softening effect results to D decreasing, 
which reveals a fact that the effective coupling J reduces 
and gives an implication with carrier-magnon decouple 
in system. 

In the conclusion, by our theoretical study the anoma- 
lous magnetization sharp drop in the vicinity of Tc comes 
from the magnon softening effect, and this softening ef- 
fect possibly comes from the carrier-magnon decoupling. 
We thanks the support of National Center of Theoretical 
Science of Taiwan. 
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